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Abstract 

We study kaon production in hot and dense hypernuclear matter with a conserved zero total 
strangeness and a conserved small negative isospin charge fraction in order to make our results 
relevant to relativistic heavy ion collisions. The baryons and kaons are treated as MIT bags in the 
context of the modified quark-meson coupling model. They interact with each other via the scalar 
mesons a, a* and the vector mesons u>, <f> as well as the iso- vector meson p. We adopt realistic sets 
of hyperon-hyperon interactions based on several versions of the Nijmegen core potential models. 
Our results indicate that the hyperons as well as the kaons are produced abundantly when the 
temperature increases and approaches the critical temperature for the phase transition to a quark- 
gluon plasma. Moreover we find that the kaons are only produced thermally and we find no kaon 
condensation in the regime explored by relativistic heavy ion collisions. 
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I. INTRODUCTION 



The quark-meson coupling (QMC) model incorporates the quark-gluon degrees of freedom 
while respecting the established model based on the hadronic degrees of freedom in nuclei. 
It describes nuclear matter as a collection of non-overlapping MIT bags interacting through 
the self-consistent exchange of scalar a, a* and vector mesons u, <fi as well as the iso- vector 
meson p in the mean field approximation with the mesonic fields directly coupled to the 
quarks In the so-called modified quark-meson coupling model (MQMC) 0, 0, 0] , it has 
been suggested that including a medium-dependent bag parameter is essential for the success 



of relativistic nuclear phenomenology. The densi 
introduced by coupling it to the scalar mesonic field 



rv-jiependence of the bag parameter is 
We have used the MQMC to study 



the properties of nuclear matter and quark deconfinement at finite temperature 0, |6|, |7| and 
the phase transition from the hadronic phase to the quark-gluon plasma j?l Q • We have also 
studied the properties of neutron stars in the context of MQMC [9]. 

In the present work we shall extend the MQMC to study the production of strange 
hadrons, especially kaons, in relativistic heavy ion collisions. Elucidating the properties of 
kaon production and kaon condensation is important to specify the tricritical point of the 
phase transition from a QGP phase to a hadronic phase. The competition between pion 
and kaon condensation in the framework of a microscopic model has been considered at 
zero and finite temperatures in the space of quark chemical potentials [10]. We consider a 
system of strange hadronic matter with conserved zero total strangeness and a conserved 
small negative isospin charge fraction in order to make our results relevant to heavy ion 
collisions. We also investigate the possibility for the onset of kaon condensation in such a 
system. 

In the mean field approximation, the scalar meson a with mass 550 MeV is supposed to 
simulate the exchange of correlated pairs of pions and may represent a very broad resonance 
observed in tctt scattering, while the vector meson u is identified with an actual meson 
having a mass of about 780 MeV. The a scalar and u vector mesonic fields couple the up 
and down flavors. Also the scalar a*(ss) and vector <p(ss) mesonic fields are actual mesons 
with masses 975 and 1020 MeV, respectively. They are introduced to couple to the strange 
flavor quark. Similarly, the mesonic iso-vector field p, which is essential for isospin charged 
nuclear matter, is an actual meson with mass 770. It exchanges the non-strange flavor 
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quarks. We assume hadronic matter to consist of the members of the SU(3) baryon octet 
and the kaon doublet. The pions are not explicitly included in the calculation because the 
contribution of one-pion exchange vanishes in spin and isospin saturated nuclear matter and 
so their inclusion will not have any effect on the quantities of interest in the present work 
such as baryon or strangeness density. The pions contribute to the total energy density and 
to the pressure, for instance, but these quantities are not calculated in the present work. 
For the purpose of calculating these latter quantities, the pions can be included as an ideal 
gas. The baryon octet is comprised of (p, n, A, E + , E°, £~, S°, H~) baryons while the kaon 
doublet is comprised of (K + , K°) mesons. In heavy ion collisions the total net strangeness of 
the system is conserved during the collision. If the quark-gluon plasma is produced during 
the collision the net strangeness for the hadronic phase may not be conserved but the total 
net strangeness for both the hadronic and quark-gluon phases will be conserved. 

Various calculations and models predict that metastable strange systems with strangeness 
fractions of order one and charge neutrality might exist in the hadronic phase at densities 
between two and three times nuclear matter saturation density 111. 12]. However, the pre- 
dicted phenomenon of metastability of strange hadronic matter and the actual values of the 
binding ener gy d epend specifically on the partly unknown hyperon potentials assumed in 
dense matter 11 , Q, ^| . Recently, unconstrained Relativistic Mean Field calculations 
of a charge-neutral strangeness-rich hadronic system have been carried out with two sets of 
hyperon-hyperon (or YY) potentials, the first set is determined from the Nijmegen hard-core 
potential Model dQ, and the second set corresponds to the potentials obtained in a re- 
cent SU(3) extension of the Nijmegen soft-core potential Model NSC9?Q H, Q, Q. The 
differences between these two sets are essentially attributable to the extremely attractive 
SS and SS interactions in the second set which allow for the possibility of deeply bound 
nuclear matter with hyperonsj^J. 

In the present work, we shall study the production of kaons in neutral as well as weakly 
isospin strange hadronic matter at finite temperature. We shall adopt realistic phenomeno- 
logical fitting parameters which are calculated from the potential models. The outline of 
the paper is as follows. In section II we write the equation of state. Then in section III we 
present the two different sets of fitting parameters for the YY potentials. In section IV, we 
present and discuss our results. Finally, we summarize our conclusions in section V. 
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II. THE EQUATION OF STATE 



The quark field i/) q (r, t) inside a bag of radius Ri representing a baryon of species i satisfies 
the Dirac equation 

i^d^ -m° q + (g qa a - g^u^) + (g q a*cr* - g q ^^) ~ 9 qP if ' t) = 0, (1) 

where m° is the current mass of a quark of flavor q. The current quark masses are taken, for 
the up and down flavor quarks, to be m u = ma = while for the strange flavor m s = 150MeV. 
Inclusion of small current quark masses for the non-strange flavors or other values for the 
strange flavor leads only to small numerical refinements of the present results. In the mean 
field approximation the meson fields are treated classically and the space like components of 
the vector fields vanish for infinite systems due to rotational invariance. As a result we get 
cj m 7 m =< ujq > 7° = c^7° and M 7 M =< O > 7° = 07° an d t • p M 7 M = r 3 < p > 7 = r 3 p7° = 
2/3^7°. The non strange (up and down) flavor quarks are coupled to the scalar er(550) and 
vector cj(780) and iso-vector p(980) mesons while the strange flavor quarks are coupled to 
(j*(975) and 0(1020). 

For a given value of the bag radius Ri and scalar fields a and a*, the quark momentum 
x qi is determined by the boundary condition of confinement which, for quarks of flavor q in 
a spherical bag, reduces to jo(x q i) = (3 q iji(x q i), where 



qi 



Vt qi {a, a*) - Rim* 



(2) 



\ Q qi (a,a*) + R,, 
We have defined the effective quark mass inside the bag as 

m *q = m q - fiV°~ - fiV'O-*, (3) 

and the effective quark energy is given by 

n qi (a, a*)/Ri = ^(xjR^ + mf. (4) 

The bag energy is given by 

E bag, = ^quarks, _ f 1 + y^ifo <T *)> ( 5 ) 



where 



jfc/ quarks i ~ ^ W 



is the total constituent quarks kinetic energy inside the bag and term is the zero-point 
energy of the quarks and Bi(a,<j*) is the medium-dependent bag parameter. 

We would like to emphasize here that we did not take into account the thermal excitation 
of the constituent quarks inside the bag 



2l|. InRefs. 



21], the total constituent 



quarks kinetic energy with thermal excitations of the populated quarks is given by 

"\ tt m (a,a*) 
q 



E. 



quarks^ 



(fq[£q-^-fli e *q + Vl. 



(7) 



where f q [e* — fi q ] and f q [e* + fi q ] are the Fermi-Dirac quark and antiquark distribution 
functions, respectively. The constituent quark effective energy is given by 



and 



Ri 



Ra 



+ (gqco^ + g q <j,<j> + -jVqpTZp) 



(fiV^ + g q<j> 4> + -j9q P T?>P) 



(8) 



(9) 



In this construction, the constituent quark chemical potential \i q is determined from the 
constraint 



n„ 



E fqK ~ Vq] ~ fq\?q + N\ 



(10) 



where n q is the number of the constituent quarks (antiquarks) inside the bag. This con- 
straint, fortunately, gives f q ~ 1 and / ~ in the actual numerical calculations. We have 
found that adopting En.© rather Eq.JTI) in the calculations doesn't affect the results given 

nnQ 

in Refs 5, 6, 7 . Therefore, the thermal excitations of the populated quarks inside the bags 

n n q 

in that construction is negligible p, 10, 17J| • 

In the simple QMC model, the bag parameter (constant) is taken as B Q corresponding to 
its value for a free baryon. The medium effects are taken into account in the MQMcP, Q| 
by coupling the bag parameter to the scalar meson fields. In the present work we use the 
following_generalized ansatz for the coupling of the bag parameter to the scalar fields a and 
a* 



Bi(a,a*) = B exp 



, / bag . bag *n 



(11) 



w ith ?^ and gf^ as additional coupling constants to be fitted from the phenomenology. 



In Ref. 



14j . the bag parameter has been considered to couple to the non strange a scalar 



field but not to the strange a* scalar field. The spurious center-of-mass energy is subtracted 
to obtain the effective baryon mass 



M * = V E ba gi 2 - < Pcm >» (12) 

where 

n q 

<Pcm>*=E4/^- ( 13 ) 
q 

The bag radii Ri and Rk for baryon i and kaon k species, respectively, are obtained through 
the minimization of the baryon mass with respect to the bag radiusjl| 

dM* 

0. (14) 



dRi 



The coupling of the scalar mean fields a and a* with quarks in the non overlapping MIT 
bag through the solution of the point-like Dirac equation should be taken into account to 
satisfy the self-consistency condition. However, this constraint is essential to obtain the 
correct solution of the scalar mean field a and a*. We think the discrepancy between our 
results in Ref. and those in Ref. jsj] is due to the coupling of the quark with the scalar 
mean fields a, a* in the framework of the point like Dirac equation. The differentiation of 
the effective hadron species mass with respect to a gives 



jp I T>a& \ _x( d <Pcm>' 
dM* ba & V &<T J a* 2 ^ da 
~da~ = M* 



(15) 



where 



da I q \ d^-qi l \ da I \ da 



d^iqi i Ri 



and 



'^bagA Arr dB t (a,a*) 

-da~ = Y 1 da (18) 
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The -S^ depends on x q {. Its solution is obtained from the point-like Dirac equation for the 
quarks and must satisfy the required boundary condition as well as the confinement on the 
bag surface jl[ It reads 

1 on n 



Ri da 



-9qc < >> (19) 



where 



*= >= n„(^-i ) + m ; fi ,/2 ' (20) 

Furthermore, the term for the spurious center-of-mass momentum correction reads 

d < vl m >i _ 1 ^ dx% 

d2 kti y ZL ) 



da R\ g da 



where 



^ = ( 2«^)+ 2i2?m^. (22) 

Similar expressions can be obtained for the differentiation of the effective hadron (baryon 
and kaon) species mass with respect to the scalar mean field a*. The zero-point en- 
ergy parameters Z\ of Eq.(4) are used to fit the actual masses of the free baryons Mi = 
939, 1116, 1193 and 1315MeV and are found to be Z i= 2.03, 1.814, 1.629 and 1.505 for 
the N, A, £ and H hyperons respectively, corresponding to a free baryon bag parameter 
Bq = (188.1) 4 MeV 4 and a free nucleon bag radius Ro = 0.6fm. The zero-point energy pa- 
rameter Zk is determined by fitting the actual mass of the free kaon = 495 MeV which 
yields the value Zk= 1.170. 

The effective Fermi-energy for baryon and antibaryon species i is given by 



<(P)= vV + M* 2 + X 4 (23) 



and 



e*W = vV + M* 2 -X 4 , (24) 



respectively, where 



Xi = [g iw u + g i4> (j> + giphip]. (25) 
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The effective Bose-energy for the kaon and antikaon species reads 



el(p) = Jp 2 + M- h : 2 + X k , (26) 



and 



r k (p) = Jp* + Mf-X k , (27) 



respectively, where 



X k = {gk^ + 9k<p<P + 9k P hkp] ■ (28) 



The indices for the kaon and antikaon species are given by k = (K + , K°) and k = (K~ , K°), 
respectively. 

The chemical potentials /ij and \x k for baryon species % and kaon species k read 

Hi = Bifi B + Si/j,s + Qi/iQ, (29) 

and 

A^fc = S k fjts + QfcyUg, (30) 

respectively. The quantum numbers Bi, Si and Qi are the baryon, strangeness and isospin 
quantum numbers, respectively. The baryon, strangeness, isospin chemical potentials /ie, 
lis and fiQ, respectively, are determined from the constraints of total baryon, strange and 
isospin densities % Eq. fpTTj) . nj°^ a ^ Eq. fjHKjl and riq ^ Eq. (f3T)|) . respectively. 
The density for each baryon and antibaryon species is given by 

dj 

Ui ~ (2tt) : 

and 

_ di 



J d 3 pf[e*(p) - f^, (31) 



d 3 P m(p) + /4 (32) 



(2tt) 3 

respectively, where di — 2 is the spin degeneracy for the octet baryon species and / is the 



Fermi-Dirac distribution function 
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Similarly, the density for each kaon and antikaon species is given by 

4 

(27T) 

and 



n» = 7 ^/<fi*K(p)-/i*] l (3L) 



^ = ?2^3 J d s pb[e* k (p)+ii k ], (35) 

respectively, where 4 = 1 is the kaon spin degeneracy, and b is the Bose-Einstein distribution 
function 

6 M = TW^TTy (36) 



When the condition for kaon condensation (KC) 22J, |23j is not satisfied, see Eq. (|52j) below, 
the kaons will not condensate and the density of the kaon condensate ti-^q = 0. In this case 
the chemical potentials ^ » s and „ Q are determined from n B , nJ M and n™*. However) 
when the condition for the kaon condensation fix = (0) is satisfied for either kaon species 
K + or K° then that species will condensate. Beyond the condensation point, the strange 
chemical potential saturates at the threshold energy for the ground states hk = €* K (0). On 
the other hand if the condition /jr = —e* K (0) is satisfied for either K~ or K° antikaon, 
then either K~ or K° antikaon will condensate. Beyond the condensation point, the strange 
chemical potential saturates at the threshold energy for the ground states hk = — e^-(O). 
We shall return to the kaon condensation when we reach Eq. (|52jl . 

The baryon density reads, 

Bs 

n B = ^2Bi[m -rii). (37) 

i 

The total strangeness density is given by 

n s = n f s + nf s } (38) 

where ng S = £)j S Si[rii — 7%] and n%f s = Y^ S Sjfcfc* — n k ] + SkUj^q- ^he total isospin 
density is given by 

n Q = n^ s + n% s , (39) 
where Uq S = Y,f S Qi[n>i ~ n%] and rig s = Y^k Qk[nk - n k ] + Qku-^q. 



The vector mean fields u, and p are determined, respectively, as follows 

Bs n Ks n 

-j [n, -«i] + L3 m ~ M + — T n KC , (40) 

lit,, lift,, lib,, 

Bs Ks 
= E "4K _ + E 2 K - n k ] + -^fn KC , (41) 

Bs Ks 

i Trip ^ m p m z 

The gw, and gt p are the meson-baryon coupling constants while g^, gk<t> and g kp are the 
meson-kaon coupling constants and are determined from the data. 

The pressure for hot and dense strange nuclear matter is equal to the negative of the 
grand thermodynamic potential density and is given by 

j Bs j „ 2 

P = 3 E j^y 3 J d3 P^ (/[<(?) - *) + MiP) + A*i]) + Pk 

1 2 2,1 2 j' 2 i 22 ^ 22 ^ 

-m.,0J H — rriiffl H — m„p w_cr i 

2" 2 2 p 2 CT 2 

where summation i runs over the 8 species of the baryon octet p, n, A, S + , S°, S" and H°, E~. 
The thermal kaon pressure is given by 

i Ks J „ l2 

^ = o E 7^3 / rf3fc / = - A**] + bfflk) + p k ]) , (44) 



i 2 2 i 2 j.2 , 22 22 -"-2*2 /<n\ 

+-m^tu z + -mffi + -m p p - -m^a A - -m a *a , (43) 



'k 2 + Ml 

where the summation k run over the kaon doublet species K + , K°. 

The scalar mean fields a and cr* are calculated by maximizing the pressure 

dP 
dP 

The pressure depends explicitly on the scalar mean fields cr, a* through the last two terms in 
Eq. (|43jl . It also depends on the baryon and kaon effective masses M* and M£ which in turn 
depend on cr and a*. If we write the pressure as a function of ({M*}, {M£}) and (a, cr*), 
the extremization of P({M*}, {M^}, a, cr*) with respect to the scalar fields a and a* can be 
written as 

dP ( dP \ dM* (dP\ 
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A similar expression can be written for J^r = 0. The coupling of the scalar mean fields a, a* 
with the quarks in the non-overlapping MIT bags through the solution of the point-like 
Dirac equation should be taken into account to satisfy the self-consistency condition. This 
constraint is essential to obtain the correct solution of the scalar mean field a, a*. 

The pressure of the kaon condensate -P^O ^ ^ ex ists, contributes to the total pressure 
of hot and dense strange hadronic matter. The total pressure of the system becomes 

ptotal = p + p KC ( 48 ) 



The contribution of the kaon condensate pressure reads 



24. 



P KC = \ UOf [4(0) - HkWM + A**], (49) 

where 

el(0)=M* k +X k (50) 

4(0) = Ml - X k . (51) 

This equation is derived by assuming that the kaon amplitude is replaced by the ansatz 
K = f '9/V2 and KK = \f 2 6 2 where / is the kaon decay constant and 6 is the dimensionless 
kaon field strength. The variation of Equation (|49|) gives 

5P KC = - A**K(0) + A**] = 0. (52) 

There are three solutions to Eq. (j52j) that determine the conditions for kaon condensation. 
The trivial solution is that kaon amplitude vanishes 9 — 0. This case means that no kaon 
condensation takes place in the system and subsequently ^KC = ®- ^ ne °ther two solutions 
are determined as follows: if [e£(0) — = either k=(K + or K°) will condensate, and if 
[ e fe(0) + fik] = either antikaon k = (K~otK°) will condensate. The density of the kaon 
condensate is given by, 

n KC = (f9) 2 (^-X k ). (53) 

In this case, the charge chemical potential becomes f^Q{e* K (0), fis) for either i^ + or K° con- 
densate or /j,q(— e* K (0), us) for either K~ or ~K° condensate since fix = Sk^s + Qk^q and 
^Tbtal j g conserV ed. It is interesting to note here that the three solutions give P^Q — 
and normally the pressure of the condensate kaons contributes indirectly via the strange 
hadronic abundances and the total energy density of the system. 
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III. FITTING PARAMETERS FOR YY POTENTIALS 



We assume that the a and u mesons couple only to the up and down quarks while a* and 
4> couple to the strange quark. We thus set g r( p = g ra * = g sa = g su = and g sp = where r 
refers to the up and down flavors while s denotes the strange flavor. By assuming the SU(6) 
symmetry of the simple quark model we have the relations g sa * = \plg ra and g S( p = —y/2g ri0 . 
The a mean field is supposed to simulate the exchange of correlated pairs of pions and may 
represent a very broad resonance observed in 7T7T scattering. We take m a = 550 MeV. The 
vector uj and p mesons are identified with the actual mesons whose masses are = 783 
MeV and m p = 770 MeV, respectively. Since the mean fields, a and u, are considered as 
< ud > condensates, they interact only with u, ci-quark in the baryons and kaons. On the 
other hand, the scalar and vector mean fields cr*,4> are considered as < ss > condensates 
and interact only with the s-quarks in the baryons and kaons. The iso-vector mean fields 
interact with the u, d-quarks in the baryons and kaons. 

The coupling of each baryon species with the vector mesons is calculated by counting the 
constituent quarks 
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9iu> = Y,9vu = 9ru, (54) 



and 



9i4> = ^29g<f> = H 9s<p, (55) 



while 



9i P = 9 gP - (56) 

With these assumptions the only free parameters left at our disposal are the quark-meson 
coupling constants g ra and g roj and the bag coupling constants g^^, gf^ for the 4 baryon 
species and these parameters are adjusted to fit nuclear properties as well as the extrapolated 
properties of hypernuclear matter. The coupling constants of the scalar and vector mesons 
to the nonstrange quarks are taken as g ra = 1 and g ruJ = 2.705 which together with a bag 
coupling constant g^^ = 6.81 yield a binding energy of 16 MeV and a compressibility Ky 1 
of 289 MeV at the normal saturation density p = 0.17fm~ 3 of nuclear matter^, ^, Q. The 
isospin vector meson coupling constant for the p-meson is taken as g qp = 8.086 to reproduce 
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the bulk symmetry energy. The coupling constants and gf^ where i = A, E, S are 

determined from the YY potentials assumptions. 

In the present calculation, we have considered two different models for the hyperon- 
hyperon interactions. Table I summarizes the values of the two sets of coupling constants 
corresponding to these two models. The basic quark-meson coupling constants as well as 
the parameters g^ 8 ^ are identical in the two models and are chosen to fit nuclear and 
hypernuclear properties. In particular, the parameters are determined by fitting the 

hyperon potentials in nuclear matter [l2j: 

U ( A N \po) = -30MeV, 
U^ N \p ) = +30MeV, 
ui N \p ) = -18MeV, 

(57) 

where the hyperon potentials are defined by 

Uf> = (M* - Mi) + {g^u + gi4> 4>). (58) 

However, the two models differ in their choices for the constants gf^- 

Model (V) is designed to mimic the consequences of the Nijmegen hard-core potential 



model D 



12 



13 



3 



It is nevertheless constrained by A and S nuclear phenomenology, and 



by a few AA hypernuclei reported to date. It accounts more realistically for the attractive 
AA and iVS interactions, but ignores the £ hyperons altogether. In this model, the medium 
constants are adjusted so that the potential of a single hyperon embedded in a bath 
of S matter becomes 

Ui(p ) = Uf(p ) = -AOMeV (59) 

in accordance with the attractive hyperon-hyperon interaction of the Nijmegen hard core 
potential Model pjl^l. Furthermore, we adopt the approximation t/#(/>o) ~ ^s(/°o) to 
fit the medium constants. The resulting U^\po/2) is about -20 MeV. 

Model (II) is designed to generate qualitatively similar baryon potentials to those obtained 
in the BHF approximation from the SU(3) extensions of the Nijmegen soft-core potential 
model NSC971S 12, 18, 3- 

The phenomenology m this model departs substantially from 
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TABLE I: Fitting parameters 

-i — i ■ , . bag bag bag bag bag bag bag bag bag 

Fit set g qa ggu g qp g N ° g AtT ta g Aa ? g Sf7 b g Sg ° g 5f7 to g Str ° g K ° g K ° g KuJ 

Model(Ia) 1.0 2.705 8.086 6.81 4.22 5.45 1.63 7.26 2.27 9.12 2.27 y/2 x 2.27 g qw 

Model(Ib) 1.0 2.705 8.086 6.81 4.22 5.45 1.63 7.26 2.27 9.12 2.27 y/2 x 2.27 l.7g qw 

Model(n) 1.0 2.705 8.086 6.81 4.22 0.0 1.63 10.28 2.27 10.17 2.27 y/2 x 2.27 g qw 

that in Model I. The NSC97 model has been tuned up to reproduce certain characteristics 
of A hypernuclei, particularly its version NSC97f. It yields particularly attractive EE, EE 
and EH interactions. Some of the shortcomings of this model are the vanishingly weak 
AA and VE interactions which are in contradiction with the little experimental evidence 
available from AA hypernuclei and from H-nucleus interaction. In this model we adjust the 
bag constants to reproduce qualitatively the binding energy curve of each hvperon 



No 



species in its own hyperonic matter B\ as those produced by the Model NSC97fQ 
binding occurs for A hyperons while E matter is deeply bound at -33 MeV per baryon at 
Pso which is twice as deep as ordinary nuclear matter, and H matter has an energy of -23 
MeV per baryon at p-o- 

The kaon coupling constants are determined from the kaon potentials when embedded 



M 122, |27J. In 



in nuclear matter = 20 MeV |2ij and = -120 to -180 MeV 

Model I, we have used two slightly different versions of these coupling constants. In set 
(la) we have taken g Kul = g qul based on simple SU (6) symmetr y w hile in set (lb) we have 
taken gxu — 1-TgVi ^ n or der to fit the kaon potential depths 26j. In both sets we take 



9k^ = 2-27 and gj^f = y/2g K ^. In Model II, we only take gx^ = g qL o based on simple 
SU(Q) symmetry. In our calculations we shall adopt gxui = gqw as the normal case and we 
shall consider the case gKw = 1-7 g qw as a special case. Here it may be worth mentioning 
that recent results |28| indicate that U^) = —50 to —60 MeV and for such a shallow kaon 
potential the renormalization factor need not be as large as 1.7. 



IV. RESULTS AND DISCUSSIONS 



We have studied hadronic matter with a conserved strangeness density ris = at finite 
temperature using the MQMC model which takes into account the quark degrees of freedom 
as well as the medium dependence of the bag constant B. We have chosen a direct coupling 
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of the bag parameter B to the scalar mean fields a and a* and the bag parameter becomes 
B = B(a,a*). The quarks and anti-quarks confined in the bags are coupled to the scalar 
mean fields a and a* and the vector mean fields u and and iso-vector mean field p. At 
first, we have studied the symmetric hot and dense strange hadronic system of baryons and 
kaons with a neutral isospin charge uq/ub = and with a conserved zero net strangeness 
ns = 0. Then, we have considered the asymmetric system with a conserved isospin charge 
nq/riB and a conserved zero net strangeness. This is relevant for heavy ion collisions where 
it is possible to have an asymmetric system with a conserved low isospin fraction in the 
range uq/ub = 0.0 to —0.2. 

The starting point of the calculation is the determination of the chemical potentials 
Pb,Ps an d Pq using the baryon, strange and isospin densities given by Eqs. ()37|) . (J3HJ) 
and (|39|). respectively. The vector mean- fields uj and <fi and the iso-vector mean field p 
are computed from Eqs.([40|). ()41|) and ([42)1 . respectively. Once the values for the chemical 
potentials (ps, Ps, Ho) are given, the values for the scalar mean fields a and a* are calculated 
by maximizing the pressure using Eqs. (j45j) and (|4"E|) . respectively. These constraints should 
take into account the coupling of the quark with the scalar mean fields in the frame of 
the points D.rao elation exact.yfll Q to obtain consistent results. The bag radius 
for each hadron species is obtained by minimizing the hadron mass with respect to its bag 
radius, Eq. (JI3jl . 

In Fig.((l}, we display the density fraction for each baryon species Xi = (n, — ni)/riB 
versus the baryonic density % at several temperatures using the parameter sets (la) and 
(lb). Furthermore we also display in the lower part of panel (a) of this figure the relative 
abundance of the kaons. The abundance of the strange hadrons up to T = 50MeV is 
completely negligible and this explains why the 5 MeV curve cannot be distinguished from 
the horizontal axes of the panels (b), (c) and (d). The strange hadrons make their first 
noticeable appearance at T « 80 MeV and they start to be produced abundantly when 
the temperature exceeds lOOMeV. The E appears to be produced more abundantly than 
the A only because it has a higher isospin degeneracy factor. On the other hand, the H 
has the lowest rate of production. The dependence of the production rates of the various 
hadronic species on the total baryon density seems to undergo a drastic change at the critical 
temperature T c = 170MeV of the QGP phase transition, becoming almost flat above this 
temperature. We also note that the modification of the kaon vector coupling which leads 
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to the two sets (la) and (lb) has very little effect on the production of the various hadronic 
species. The effect can be mainly seen above T c where set (la) produces slightly more H and, 
to a lesser extent, more £ and kaons than set (lb) because of the stronger vector coupling 
for kaons in (lb). Kaon production in the system increases with temperature especially for 
temperatures exceeding T « 100 MeV. We see that a large mesonic strangeness n$ > 0.2 can 
be reached at the critical temperature T c = 170 MeV when kaons are produced abundantly 
in the system. 

In Fig.©, we compare the baryon and kaon abundances obtained using parameter sets 
(la) and (II). The phenomenology in model (II) departs substantially from that in (la) for 



large baryonic and strangeness densities where the phase transition to EE takes place[ 
We see that at high temperatures the A is more abundantly produced with the set (la) 
than when the set (II) is used. The X has the opposite behaviour, as its abundance ob- 
tained from the set (la) is less than that from the set (II). At the same time, there is no 
noticeable difference between the two models in the nucleon, kaon and cascade production 
rates. Furthermore, there is not much difference in total strange baryon production. Gen- 
erally speaking, it can be concluded that kaon production and condensation in the context 
of MQMC is not likely to be affected drastically by which set is being used from amongst 
the different sets of parameters extracted from phenomenology. 

In Fig.Q, we display the anti-baryon to baryon abundance ratio for each baryon species 
Hi/ rii at several temperatures as obtained by using parameter set (la). It is seen that 
for temperatures T < lOOMeV the production of the anti-baryons is negligible. However, 
for higher temperatures there is significant anti-baryon production particularly when the 
temperature reaches and exceeds T c . Furthermore, we notice that the anti-baryon abundance 
ratio T%i/ Hi, is particularly large at small baryon density This ratio also seems to increase 
with the baryon strangeness number as follows: 

S AS AT 

I > A-E > F (60) 
This enhancement of anti-hyperon abundance ratio over the antiproton abundance ratio in 
the relativistic nuclear collisions has been observed earlier j^. Also included in Fig.Q are 
the results obtained in the ideal gas approximation which assumes noninteracting hadrons 
in a thermal bath. In this approximation, the hadrons do not modify their masses, do not 
interact with each other and there are no mean fields at all. Although this approximation 
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is too simplistic, it can be useful to obtain rough first estimates. It is seen from Fig.® 
that the ideal gas results greatly underestimate the anti-baryon production, especially at 
high temperature or at high density. In Fig.(J3J), we show the anti-kaon to kaon abundance 
ratio versus the baryon density at various temperatures using parameter set (la) and 
in the ideal gas approximation. We note that the anti-kaon abundance ratio increases with 
temperature. We also note that again the ideal gas approximation gives much lower estimates 
for the anti-kaon to kaon abundance ratio than the MQMC, especially at high temperature 
or at high density. 

Fig. © examines the crucial issue of the onset of kaon condensation in heavy ion collisions. 
Kaon condensation implies that the equation of state will be softer and that more hyperons 
and kaons will accumulate in the system. Indeed, it is known that kaon condensation plays 
a significant role in the physics of neutron stars |27(. One of the major themes of the present 
work is to study the possibility for the onset of kaon condensation in heavy ion collisions. 
The necessary condition for the onset of kaon condensation is to nontrivially satisfy Eq. (|52|) . 
For the K° or K + kaons, condensation is possible only when 

Hk = 4(0) (61) 
while for the the antikaons TC° or K~, condensation is possible only when 

Hk = -4(0). (62) 

When the condition for the onset of kaon condensation is satisfied, the kaons start to con- 
dense smoothly and the strange chemical potential could be evaluated either from Eq. (j61|) 
or (|H2*|) . After the kaon condensation starts to take place in the system, the strange chemical 
potential fis becomes an input parameter and the amplitude of the kaon condensation is 
calculated by conserving the total strangeness of the system. The contribution of the kaon 
condensate to the total density of the system is calculated from Eq.()53J). To examine the 
possibility for the onset of kaon condensation in our model we plot in Fig.© the kaon chemi- 
cal potential \ik and the kaon threshold effective energy e£(0) versus the baryonic density % 
at various temperatures. The isospin chemical potential fiQ vanishes for symmetric strange 
hadronic matter so that in this case Eq. (|HUj) reduces to fix = Us where /is is the strangeness 
chemical potential. We see from panel (a) of Fig.© that the kaon chemical potential hk 
is always less than the kaon threshold energy e£(0) for the all baryon densities n B so that 
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Eq.(|61|) can never be satisfied. A similar conclusion can be drawn for the anti-kaons and 
Eq.([62j) as can be seen from panel (b) of the same figure. This indicates that kaon conden- 
sation does not take place at any temperature and density in heavy ion collisions where the 
total strangeness is conserved to zero for a system that consists of baryons and kaons. The 
kaons are only produced thermally in this system. The situation is, of course, completely 

different in neutron stars where there is enough time for the conservation of strangeness to 

i ■ " 

be violated by the weak interaction [2jj. 

In Figs.© and (J2J), we display the baryon and kaon effective masses versus the baryon 
density ng. We note that, in the context of the MQMC model, the effective masses do 
not vanish at large baryon densities, unlike the standard Walecka model where the effective 
hadronic masses tend to vanish at large baryon density. The realistic case of effectively 
massive hadrons at all densities is actually a characteristic result of the MQMC. Indeed, in 
the hadronic phase, the hadron's effective mass should not vanish until chiral symmetry is 
restored or a QGP phase transition takes place. At the phase transition, there is a sharp 
drop in the hadron's effective masses and all hadrons in the system will dissolve into their 
constituent quarks. Furthermore, the quarks attain their current masses which are very 
small for the up and down flavors. At this point, however, bubbles of the QGP are formed 
in the system. Therefore, we do not find it is appropriate, in the realistic situation relevant 
to heavy ion collisions, to extrapolate the study of the onset of kaon condensation to massless 
hadrons since these massless hadrons can exist only at the point of the phase transition. 

We have also extended MQMC to study hot and dense asymmetric strange hadronic 
matter with a conserved total zero strangeness ris/ris = and a conserved finite negative 
isospin density fraction uq/ub < 0. Generally speaking, the negative value of the isospin 
density fraction uq/ub reflects the fact that the number of neutrons is more than the number 
of protons (i.e. n n > n p ) in the initial conditions of the colliding heavy ions. The results for 
the asymmetric system are therefore more relevant to relativistic heavy ion collisions. We 
present our results for such a system using the parameter set (la) in Figs. (JBJ) and (jHJ) for 
hq/tib = —0.10, while the results for uq/hb = —0.20 are shown in Figs.^HJ) and (fTTf . 

In Figs.© and (JTUJ), we display the net kaon abundance ratio xk = {nK—n K )/n B and the 
ratio of anti-kaon abundance over kaon abundance for the charged and neutral kaon species 
(K + ,K°) at several temperatures for the isospin charge density fractions uq/hb = —0.1 
and uq/ub = —0.2, respectively. In panel (a) of both figures, we display the net kaon 



18 



abundance versus the baryon density %. It is seen that xko = {n^o — n-^o )/ns is much 
larger than %k+ = (n K + —n K -)/riB in such a system, especially at higher temperatures. The 
conservation of net strangeness for the asymmetric system is more easily preserved by neutral 
kaons at high temperatures. Consequently, when the negative isospin uq/ub is increased, 
the neutral kaons become more dominant than the charged ones. ^From Fig.fjSJb), we see 
that the K° and K + are produced more abundantly than their anti-particles for the system 
with lower asymmetry. However, when the asymmetry is increased, we see from Fig. fllUl b) a 
dramatic change for the charged kaons with the anti-particle K~ becoming more abundant 
than the K + at high densities in order to conserve isospin. Actually, in addition to the kaons, 
we have found that the baryon species E~, S°, S + also play an important role in conserving 
the isospin uq/ub at high temperature when the isospin charge fraction is increased. The 
proton density fraction tends to saturate at some value when the negative isospin charge 
fraction becomes rather large (uq/ub = —0.2,n p /nB ~ 0.30) while the density fraction 
for neutrons continues to increase with temperature. However, at higher temperatures, 
the conservation of isospin charge fraction is preserved by creating more charged strange 
hadrons. 

We have also examined the possibility for the onset of kaon condensation in the asym- 
metric system. The necessary condition for kaon condensation is that the kaon threshold 
energy e£(0) equal its chemical potential \ik- In FigsQ and (JTTJ), we display the kaon chem- 
ical potential and the kaon threshold effective energy e£(0) versus the baryonic density 
tib at various temperatures and for uq/ub = —0.1 and uq/ub = —0.2, respectively. We 
see that the kaon's threshold energy is always larger than its chemical potential, and their 
curves become almost parallel at large baryonic density. Furthermore, when the nuclear 
matter asymmetry is increased by increasing the negative isospin uq/ub to reach the maxi- 
mal value of uq/ub = —0.5, appropriate for neutron matter, the onset of kaon condensation 
is not possible even at very large baryonic densities and we find that the strange chemical 
potential us in the medium is always modified in such a manner as to avoid the onset of kaon 
condensation. These results show that no signature can be expected for kaon condensation 
in heavy ion collisions and kaons can only be produced thermally. 

We have also examined how the production of the charged and neutral kaons and their 
anti-particles at high temperature varies with respect to the isospin charge ratio uq/ub- 
In fig.(|12|). we display the kaon and antikaon abundances versus the baryonic density % 
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at temperature T = 170 MeV for uq/ub = -0.1, -0.2, -0.3. The temperature T = 170 
MeV is particularly interesting because it is accessible at RHIC energies and is relevant 
to the freezeout chemical potential 0, [nj]. Part (a) of fig.()12jl displays the charged and 
neutral kaon abundant density fractions, xk, while part (b) displays the densities for charged 
(upper panel) and neutral (lower panel) kaons and their antiparticles. It is seen from part 
(b) that the kaons and their antiparticles are produced abundantly at high temperatures 
even at very low baryon densities. However, the net density fraction for charged kaons 
Xk = (n>K+ ~ n-K-)/n>B tends to be small and varies weakly with the baryon density 
Furthermore, it decreases when the asymmetry of the system increases by increasing the 
negative value of the isospin charge fraction uq/hb- This reflects the fact that the negative 
charged kaons become more dominant as the magnitude of the isospin charge asymmetry 
increases and exceeds Iuq/ubI = 0.25. Here it must be stressed that in a realistic situation 
the isospin charge fraction of the normally used colliding heavy ions lies in the range uqIub = 
to —0.1 and consequently K + production is expected to be more dominant than K~ in 
RHIC. However, more asymmetric nuclear matter may be reached experimentally in the 
future by using different isotopes for the colliding heavy ions. Generally speaking, the net 
production of charged kaons is relatively smaller than that of the neutral kaons, indicating 
that the neutral kaons play a more significant role in asymmetric nuclear matter in order to 
conserve the net zero strangeness of the system. 

Fig. ()13|) displays the variation of kaon and antikaon abundances with respect to the tem- 
perature for several values of the isospin charge fraction uq/ub- The calculations are carried 
out at the baryonic density ps = 0.05fm~ 3 . This rather small baryonic density is accessible 
in RHIC and has a special interest at high temperature in particular in studying the freeze- 



out chemical potential 



3d, El- The phase transition from the hadronic phase to QGP is 



expected to take place at this low baryon density and temperature T = 170MeV. It is seen 
that kaons are produced abundantly when the temperature exceeds T > 120MeV and this 
production becomes more significant when the temperature reaches T = 170 MeV. Further- 
more, it is seen that the neutral kaon density, %o, lies always above the antikaon density 
n-^o but this is not the case for the charged kaons, and tlk- ■ The net density fraction 
for the charged kaons decreases as the magnitude of the negative isospin charge asymmetry 
increases while net density fraction for the neutral kaons has the opposite behavior. More- 
over, the net density for the neutral kaons is always greater than that for the charged kaons 
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for asymmetric nuclear matter. This reflects the fact that the K° abundance is dominant at 
all temperatures and more so as the magnitude of the isospin charge asymmetry increases. 

In Fig.flHJ), we display the baryon, strangeness and isospin chemical potentials versus 
the baryonic density % at T = 170MeV for the symmetrical and asymmetrical nuclear 
matter. We also indicate on this figure the estimated freezeout chemical potential at SPS 
and RHIC energies. In particular it is known that freezeout in RHIC occurs at the small 
chemical potential fis ~ 46 ± 5MeV and at a temperature T « 174 ± 7MeV (30! . The 
substantial decrease of the baryon chemical potential from \xb ~ 270MeV at SPS energies to 
fiB ~ 45MeV at RHIC energies shows that, at mid-rapidity, we are dealing with a low net 
baryon density in the medium [30!] . In our model, the baryon chemical potential //# ~ 50MeV 
for the critical temperature T « 170MeV corresponds to the very small baryonic density 
riB ~ 0.015fm~ 3 . Other models, however, may lead to somewhat higher values for the 



baryonic density for the chemical freezeout potential 



and we shall refer to a baryon 



density % ^ 0.05fm 3 when we speak about RHIC physics. 

V. CONCLUSIONS AND SUMMARY 

We have studied strange hadronic matter composed of the baryon octet and kaon doublet 
with a conserved zero total strangeness us/ub = ^n^ ar ^ 0nS + ^Kaons^ _ g anc j a 
conserved finite isospin fraction uqIub- RHIC physics is reached at the rather low baryonic 
density % ~ 0.05fm -3 and high temperature T rs 150 — 170 MevjscJ. We have not found 
any evidence for kaon condensation even in highly asymmetric dense nuclear matter. The 
kaon threshold energy is always higher than the kaon chemical potential given by Eq. (J30|) 
so that the kaons are only produced thermally in the system. The production of the strange 
hadrons starts smoothly at low temperatures and then they are abundantly produced when 
the temperature reaches T « 170MeV. 

We have solved the baryon, strange and isospin chemical potentials UbiUs an d fiQ, re- 
spectively, self-consistently. In the hot or even warm dense medium, the strange chemical 
potential ^5 is finite (nonzero) even for a system with zero total strangeness and it likely 
modifies itself in a manner to decrease the kaon chemical potential to be always lower than 
the kaon threshold energy. The results of other microscopic models also disfavour kaon with 
respect to pion condensation jloT] . On the other hand, in cold dense matter, the strange 
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chemical potential /j,s is normally set to zero for a system with zero total strangeness. In 
this case, it is possible for the antikaon (kaon) to condense at zero temperature at very high 
baryonic density. Therefore, the results at finite temperature suggest that /is must be actu- 
ally assigned a small finite value at zero temperature in order to maintain the conservation 
of the zero total strangeness for cold dense matter. Such a modification would also reduce 
the possibility for the onset of kaon condensation at T = 0. 

It is interesting to note here that the effective hadronic masses do not vanish in the 
hadronic phase in our model. The effective hadronic masses attain small values only at 
the chiral phase transition and the phase transition to the QGP. In the QGP, the hadrons 
dissolve and disappear and the kaons will not be produced. On the other hand, at very 
high densities and low temperatures the quark color superconductivity will dominate. The 
strange hadrons and their anti-particles are normally produced at high temperatures and 
low baryonic densities which are the ideal circumstances at RHIC energies. Furthermore, 
when the system freezes out, the strange baryons and kaons might survive with the baryonic 
sector having a negative net strangeness n^ aryOI1S /n£ = 0.0 to —0.25 and the mesonic 
sector (i.e. kaons)having a positive net strangeness, but the total strangeness for the baryons 
and kaons is conserved to zero. Moreover, although K production always dominates over 
K , the K~ becomes more abundant than K + for highly asymmetric nuclear matter (i.e.: 
"^q/ub < — 0.2), especially at high temperatures or large baryonic densities. 
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FIG. 1: The density fraction for hadron and anti-hadron species x« = (rij — n^jriB versus baryonic 
density at several temperatures using the parameter sets (la) and (lb), a) Nucleons and Kaons, 
b) A, c) E, d) S. 
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FIG. 2: Same as Fig.Q but using parameter sets (la) and (II). 
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FIG. 3: The ratio of anti-baryon to baryon abundance ni/rii for each baryon species versus 
baryonic density at several temperatures using the parameter set (la) and for the ideal gas 
approximation, a) N, b) A, c) S, d) S. 
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FIG. 4: The ratio of anti-kaon to kaon abundance n^/n^ at various temperatures using the 
parameter set (la) and for the ideal gas approximation. 
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FIG. 5: The strange chemical potential fj,s and the kaon threshold effective energy e£(0) versus 
baryonic density ns for various temperatures using the parameter set (la), a) Kaons, b) Anti-kaons. 
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FIG. 6: The effective mass M* for each baryon species versus baryonic density ng for various 
temperatures using the parameter set (la), a) N, b) A, c) S, d) S. 
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FIG. 7: The effective kaon mass M£- versus baryonic density n# for various temperatures using 
the parameter set (la). 
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FIG. 8: The kaon abundance versus baryonic density ub for various temperatures for uq/ub = 
—0.1 using the parameter set (la), a) The kaon density fraction xk = (jik — nK)/nB, b) The ratio 
of anti-kaon abundance over kaon abundance nx/nx- 
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FIG. 9: The strange chemical potential fj,s and the kaon threshold effective energy ejjl(0) versus 
baryonic density at various temperatures for uq/ub = —0.1 using the parameter set (la), a) 
Kaons, b) Anti-kaons. 




FIG. 10: Same as Fig.® but with n Q /n B = -0.2. 
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FIG. 11: Same as Fig.® but with n Q /n B = -0.2. 
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FIG. 12: The kaon and anti-kaon abundances versus baryonic density ub for various isospin charge 
ratios hq/ub at temperature T = 170 MeV using the parameter set (la), a) The density fraction 
xk = (wk — n-j^j/nB for the charged and neutral kaons, b) The kaon and anti-kaon abundances 
for charged (upper panel) and neutral (lower panel) kaons. 




FIG. 13: Same as Fig. <jl2j) but versus temperature at the baryonic density ns = 0.05fm 
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FIG. 14: The dependence of baryon, strange and isospin chemical potentials hb, MS an d MQj 
respectively, on the baryon density ns at T = 170MeV for a symmetric and an asymmetric matter. 
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